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ABSTRACT 

We have performed a 1.3 cm survey of 24 massive young stellar objects 
(MYSOs) using the Expanded Very Large Array (EVLA). The sources in the 
sample exhibit a broad range of massive star formation signposts including In- 
frared Dark Clouds (IRDCs), UCHII regions, and extended 4.5/um emission in 
the form of Extended Green Objects (EGOs). In this work, we present results 
for G35.03+0.35 which exhibits all of these phenomena. We simultaneously im- 
age the 1.3 cm NH3 (1,1) through (6,6) inversion lines, four CH3OH transitions, 
two H recombination lines, plus continuum at 0.05 pc resolution. We find three 
areas of thermal NH 3 emission, two within the EGO (designated the NE and 
SW cores) and one toward an adjacent IRDC. The NE core contains an UCHn 
region (CM1) and a candidate HCHn region (CM2). A region of non-thermal, 
likely masing NH3 (3,3) and (6,6) emission is coincident with an arc of 44 GHz 
CH3OH masers. We also detect two new 25 GHz Class I CH 3 OH masers. A 
complementary Submillimeter Array 1.3 mm continuum image shows that the 
distribution of dust emission is similar to the lower-lying NH3 lines, all peaking 
to the NW of CM2, indicating the likely presence of an additional MYSO in this 
protocluster. By modeling the NH 3 and 1.3 mm continuum data, we obtain gas 
temperatures of 20-220 K and masses of 20-130 M . The diversity of contin- 
uum emission properties and gas temperatures suggest that objects in a range of 
evolutionary states exist concurrently in this protocluster. 
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Subject headings: stars: formation — stars: massive — ISM: individual objects 
(G35.03+0.35) — techniques: interferometric 



1. Introduction 

Massive star formation is a phenomenon of fundamental importance in astrophysics, 
yet it remains poorly understood. Massive protostars form in complex clusters and predomi- 



natel y at distances greater than a kiloparsec making them challenging to study (IZinnecker and Yorke 



20071 ). Unfortunately, the tools used to study low mass YSOs (primarily near-IR imaging) 
are largely inapplicable due to extreme dust obscuration even into the mid-IR. Radio and mil- 
limeter wavelengths penetrate the dust, and have revealed a wide variety of phenomena asso- 
ciated with MYSOs: CH 3 OH, H 2 0, and OH masers, hyper-compact (HC) and ultra-compact 
(UC) Hll regions, recombination lines, infrared dark clouds (IRDCs), warm (> 30 K) dust 
cores, massive outflows, extended 4.5 //m emission from shocks, and hot core line emission. 
However, these signposts have been compiled from a heterogeneous set of observations with 
varying angular resolution and sensitivity, making correlation analyses difficult. 

In order to advance our understanding of MYSOs, we are using the EVLA to observe 
24 MYSOs with ~ 10, 000 AU resolution in 1.3 cm continuum and a comprehensive set of 
diagnostic lines simultaneously. This resident shared risk observing project is among the 
first to utilize 16 spectral windows with high spectral resolution including: the NH 3 ladder 
from (1,1) to (6,6) (Ei=23 - 540 K); four CH3OH transitions that can show maser emission; 
two H recombination lines; and other potential hot core species; p lus reasonable cont i nuum 



sensitivity. The majority (22/24) of our sample originates from the ICyganowski et al.l (120081 ) 
catalog of MYSO outflow candidates (extended green objects, EGOs, selected based on 
ext ended 4.5/mi emission) . Many EGOs are located in IRDCs (for an overview of IRDCs 



sec 



Rathborne et al.l 120061 ) . with a smaller subset associated with HC or UCHII regions. 
To explore the range in the properties of MYSOs associated with extended 4.5/xm shock 
emission, the sample includes five EGOs with compact HII regions, as well as two non- 
EGOs for comparison (an IRDC and a UCHII region). In this paper, we present the first 
results from the survey by describing the observations of one source in detail: G35.03+0.35 
(hereafter, G35.03). To show the location of dust emission with respect to the centimeter 
emission, we include complementary 1.3 mm Submillimeter Array (SMA) observations. 
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1.1. Background on G35.03+0.35 



G35.03 is an EGO at a distance of ~ 3.4 kpc with a bipolar 4.5 /iin morphology ori- 
ented NE-SW at a systemic v elocity of ~ 53.1 km s _1 based on H 13 CO + (3-2) single dish 
data (jCyganowski et al.1 120091 ) . A saturated Spitzer MIPSGAL 24 /mi source is located 
near the center of the two lobes and is coincident wit h five compact centimet er continuum 
sources called CM1...CM5 in order of decreasing flux ( jCyganowski et al.ll201ll ). CM1 has a 
flat centimeter spectral index (a ~ — 0.1, S v oc v a ) consistent with optically thin free-free 
emission from an UC Hll region (also see lKurtz et al.lll994t lArgon et al.l 120001 ). In contrast 
CM2 h as a rising spectrum {a ~ +0.7), suggesting it may be a HC Hll region or wind/jet 
source (jCyganowski et al.l 120111 ). Coinci dent with CM2 are blue shifted Class-II 6.7 GHz 
CH3OH masers (Ay ~ 6 - 12 km s" 1 ; ICveanowski et al.ll2009f ). blueshifted OH masers 
(AV ~ 3 - 14 km s" 1 : lArgon et al.ll2000h . and redshifted H 2 masers (Ay ~ 13 — 17 km s 1 ; 



Forster and Caswelllll999l ). Additionally, an arc-like structure of outflow -tracing Class-I 44 



GHz CH3OH masers are located toward the Eastern side of the EGO ( jCyganowski et al. 
20091 ). 



Observations 



The NRAcB Expanded V ery Large Array (EVLA) observations are summarized in Ta- 



ble [U see iPerley et al.l ( 1201 if ) for additional details on the EVLA system. The EVLA data 



were calibrated and imaged in CASA@. All images were restored with a synthesized beam of 
3'!7 x 3f0 (P.A.=-50°), corresponding to ~ 12,600 x 10,000 AU at 3.4 kpc. 

The SMA|f| observations are also summarized in Table 1. The sidebands were centered 
at 220.1 and 230.1 GHz. The data were calibrated in MIRIAD, then exported to CASA 
for self- calibration and imaging. In this paper we only consider the 1.3 mm continuum 
data estimated from line-free channels in the uv-plane. The SMA spectral line data will be 
described in a future publication. 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under agreement by the Associated Universities, Inc. 

2 http://casa. nrao.edu/ 

3 The Submillimeter Array (SMA) is a collaborative project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy & Astrophysics of Taiwan. 
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3. Results 

As shown in Figure [lb,, we detect three major regions of NH3 (1,1) emission: (i) the 
strongest emission is coincident with the central region of the bipolar 4.5 fxm nebulosity, 
the five centimeter continuum sources described in §1.1, and 6.7 GHz CH 3 OH masers; (ii) a 
slightly weaker and more compact core coincident with the SW lobe of the bipolar 4.5 jj,m 
nebulosity; and (iii) a weak but extended core coincident with the IRDC ~ 30" SW of the 
EGO. Figures [Tb and c show the intensity weighted velocity field and velocity dispersion of 
the NH3 (1,1) emission, respectively. The emission spans a velocity range of about 4 km s _1 , 
with a distinct NE-SW blue-red velocity gradient across the EGO. This velocity gradient is 
not clearly associated with an outflow, but instead seems to trace a shift in systemic velocity 
across the region. The mean systemic velocity toward the EGO and the "Dark" core to 
the SW is ~ 54 km s _1 . In the vicinity of the "NE" core (and the centimeter sources), the 
velocity dispersion shows a marked increase from ~ 0.6 km s _1 to values as high as ~ 1.9 
km s _1 . There are also "fingers" of weak NH 3 emission evident around the periphery of the 
EGO - these fingers are not apparent in the integrated intensity map (Fig. dh) because as 
demonstrated in Fig. [lb they have extremely narrow velocity dispersion (~ 0.25 km s _1 ). 

Figure [2b shows a detailed comparison of the EV LA NH3 (1,1) integrated intensity, SMA 



1.3 mm continuum, and the VLA 3.6 cm continuum (jCyganowski et al.ll2011l . CM1 and CM2 
are also detected in the line free channels of the current EVLA data; albeit at lower angular 
resolution and sensitivity). Overall, there is strong agreement between the NH3 and 1.3 mm 
continuum emission. The strongest 1.3 mm emission arises from the NE core, coincident 
with CM1 and CM2. Though it is difficult to quantify the free-free vs. dust contributions 
at 1.3 mm due to the mismatch in the continuum resolution (and uv— coverage) between 
the SMA, VLA, and EVLA (not shown), it is clear that a significant fraction must be due 
to dust. The morphology of the NH 3 (1,1) and (2,2) emission (see Fig. [2b), along with the 
1.3 mm emission to the NW of CM2, suggests the presence of a third compact source in the 
NE core lacking in free-free emission. Millimeter continuum emission is also visible to the 
East-SE of the NE-core that partially coincides with the arc of 44 GHz CH 3 OH masers, and 
mimics the morphology of the NH3 (1,1) integrated intensity. The SW NH3 core also has 
a 1.3 mm counterpart. The "Dark" ammonia core was not detected by the SMA, probably 
due to the smaller primary beam of these data. 

Figure [2b gives a detailed comparison of the NH 3 (1,1), (3,3), and (6,6) integrated in- 
tensities. Significant (3,3) emission is detected toward both the NE and SW core regions, 
though the strongest emission is compact and located East of the EGO, toward the SE termi- 
nus of the 44 GHz maser arc. Interestingly, this region shows only weak para-NH3 emission 
suggesting a non-thermal origin for the (3,3) emission, though it appears to be connected 
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to CM2 through a weaker ridge of warm thermal NH3 emission. The (6,6) transition is also 
detected at the (3,3) peak, as well as toward the candidate HC Hll region CM2. 

We also detect strong, spatially-unresolved emission in the CH 3 OH 5 2 — 5i and 8 2 — 81 
lines at two distinct positions (Fig. [3]). Using an upper limit to the fitted size of the emission 
(F/5 x0'.'6), their peak flux densities (187±2 mJy beam -1 and 39±2 mJy beam -1 ) correspond 
to brightness temperatures (T&) exceeding 400 and 85 K for the stronger and weaker source, 
respectively. These high T& compared to the energy levels above ground (57 and 105 K) 
strongly suggest that these transitions a re masing (25 GHz CH3OH masers have also been 



observed for example in IRAS16547-4247. IVoronkov et al.ll2006l ). Interestingly, the transition 



of peak maser intensity switches from 5 2 — 5i for the stronger maser spot, to 8 2 — 81 for the 
weaker maser spot. The positions and peak velocities are (J2000) 18 h 54 m 01.043 s (±0.002s), 
+02°01'16 / .'86 (±0"02) at 52.8 km s -1 , and 18 h 54 m 00.562 s (±0.002s), +02°01'17"67 (±0"02) 
at 55.2 km s _1 for the strong and weak masers, respectively (see Fig. [2^,b,c). 

No emission greater than 5a was detected for the S0 2 , DC3N, HC5N, and HC9N tran- 
sitions. The H63a and H64a recombination lines were detected toward the CM1 UC Hll 



region at a peak velocity of 55.9 km s -1 . Following the method of iGaray et al.l (119861 ). we 
use the line-to-continuum intensity ratio of 0.46, the fitted FWHM linewidth of 17.6 km s _1 , 
and the diameter of CM1 determined from the EVLA 1.3 cm continuum (1'.'3) to derive an 
electron temperature of 7900 K and density of 1.3 x 10 4 cm -3 . 



3.1. NH 3 Temperatures 



Over the last few decades considerable progress has been made in understanding th e 
excitation of the NH3 molecule under astrophysical conditions (see e.g. iMaret et al.l 120091 ). 
For example in cold regions, analysis of the (1,1) and (2,2) transitions has evolved toward 
simultaneous least squares-fitting of the line profiles - freeing the analysis fro m previous 



failures wherever the intensity of (2,2) rivals (1,1) (e.g. iRosolowsky et al.l 120081 ). However, 
this type of analysis typically assumes that the kinetic temperature is much less than the 
energy gap between levels (41.5 K) such that only the (1,1) and (2,2) rotational levels are 
populated. Due to these assumptions, this technique becomes less accurate above tempera- 
tures of about 30 K. Generic molecular modeling packages exist that can be used to model 
NH 3 emission from higher lying transitions using either radiative transfer with collisional 
excitation, or assuming LTE. However, no package currently combines separate ortho and 
para collision rates (with realistic collision partners) with the detailed hyperfine structure of 
NH3 and a partition function that includes all relevant states that may be excited in warm 
gas (though we can expect this situation to improve in the near future). 
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For the current preliminary analysis we have pursued the following course: (1) in 
the colder regions of G35.03 we have fit simultaneously the full hyperfine structure of the 
(1,1) and (2,2) lines using a non-linear, least squares Gaussian fitting routine given a sin- 
gle common line-of sight- velocity, line width, excitation, and kinetic temperatures in IDL 



( iFriesen et al.l 120091 ); (2) In the warmer, and coincidentally more kinematically complex re- 
gions we have used the LTE method in the CASSIS0 package to model all six observed NH 3 
transitions, including multiple velocity components where necessary To compare methods, 
we also fit the cooler regions with CASSIS. Figure [3^ shows the results from technique 
(1), with grey areas indicating where this method's assumptions are invalid. The fitted pa- 
rameters for both techniques toward several representative positions are given in Table |2j 
The corresponding spectral line fits for these positions using technique (2) are shown in 
Figure [HH-h. 

As shown in Table |2j the temperatures in the kinematically-simple Dark and SW core 
regions are in the 20-30 K range, with the latter being somewhat warmer on average. Both 
methods do a reasonable job here, with method (2) predicting about twice the total column 
density of (1), likely due to its more complete partition function. Toward the center of the 
NE core and along the 44 GHz maser arc, technique (1) breaks down for two reasons: (i) 
the gas is warm (> 30 K) throughout this region (or even non-thermal; Fig. [3b) and (ii) 
the lines become more kinematically complex with at least two distinct velocity components 
separated by about 1.7 km s _1 , or in the case of CM2 what appears to be a significantly 
blueshifted hot (220 K) outflow component. Towards CM1, the redder velocity component 
shows temperatures as high as ~ 70 K. 



3.2. Dust Mass 



The average c olumn densitie s and m asses of the NE and SW cores based on the method- 
ology described in iBrogan et al.l (120091 in Eqs. 4 and 5) using the 1.3 mm SMA data are 
shown in Table [2J The integrated 1.3 mm flux densities of the NE and SW cores are 831 
and 170 mJy with estimated sizes of 11.0" and 6.8", respectively. For the NE core, we have 
subtracted 25 mJy to account for the 3.6 cm flux of CM1 and CM2 (see Figj2]) assuming 
they are d ue entirely to free- f ree em ission and using the 1.3 cm to 3.6 cm spectral indices 
derived by ICyganowski et al.l ( 120111 ). For the range of temperatures derived from the NH 3 
fits described in §3.1, the NE and SW cores have dust-based gas masses in the 50 - 132 M 
and 22 - 45 M range, respectively, and average column densities of (3 — 7) x 10 



22 —2 

~ z cm . 



4 CASSIS has been developed by CESR-UPS/CNRS ([http://cassis.cesr.frl 
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For comparison, iHill et al.l ( 120051 ) report a SEST SIMBA 1.2 mm based mass for the whole 
G35.03 region of 390 M assuming a dust temperature of 20 K and a total size of 55" (in- 
cludes the "Dark core" in addition to the NE and SW cores). The Bolocam Galactic Plane 
Survey (BGPS) catalog report s a 1.1 mm flux of 1380 and 3030 mJy for G35.03 within 40" 
and 80" apertures, respectively ( iRosolowsky et al.ll2010l ). Assuming Td us t = 30 K for r < 20" 
and T dust = 20 K for 20" < r < 40", we find a total mass for the region of 350 M Q (550 
if we multip ly by the BGPS "correction factor" of 1.5), in reasonable agreement with 



Hill et al.l ( 120051 ) and our SMA results, accounting for both the smaller primary beam and 



spatial filtering of the interferometer. 



4. Discussion 



4.1. Protocluster Nature of G35.03+0.35 



Our high angular resolution observations reveal the presence of at least three MYSOs 
(possibly four if CM3 is included) within 20,000 AU of one another in the NE core (CM1, 



Hunter et al. 


2006; 


Rodon et al. 


2008) 



Component CM1 is a modest UC H II region associated with warm NH3 and hosts two ve- 
locity components which may indicate further unresolved structure. CM2 appears to be a 
HCHll region or wind/jet source exciting OH, H 2 0, and Class II CH 3 OH masers, and a hot 
(220 K) blueshifted outflow. The millimeter continuum source NW of CM2 exhibits the peak 
dust emission, no centimeter continuum emission, but the broadest Gaussian (non-outflow) 
NH 3 profiles. Such diversity among the continuum properties and in the molecular gas 
temperatu res toward these different objects is an increasingly familiar patter n seen in pro 



tochi sters (IZhang et al.ll2007t Brogan et al.l 120071 ). Similar to NGC6334I(N) (IBrogan et al. 
20091 ) . the systemic velocities of these objects differ by up to 1.7 km s _1 , providing a measure 
of the cluster dynamics. However, detailed study of the kinematics requires higher angu- 
lar resolution followup observations, in particular to spatially resolve the multiple velocity 
components toward the NE core objects. 



5 BGPS: 



http: / /irsa.ipac.caltech.edu/data/BOLOC AM.GPS / 
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4.2. New Masers in G35. 03+0.35 



As Figs. [2b, and|3]3 demonstrate, the (3,3) and (6,6) NH3 emission toward the terminus 
of the 44 GHz CH3OH maser arc shows significant departure from LTE. Unfortunately, our 
current angular resolution precludes the measuremen t of the high brightness tempera tures 
required for absolute confirmation of maser emission. IWafmsley and Ungerechtsl ( 119831 ) first 
predicted that the (3,3) transition could undergo weak population inversion for densities 
of ri ~ 10 4-5 cm -3 and temperatures of ~ 50 K - physical conditions common in massive 
molecular outflows an d similar to the conditions required to excite 44 GHz CH3OH masers 
( jVoronkov et al.l 120051 ). We speculate that the non-thermal NH 3 emission traces a bow 
shock from an outflow originating from CM2. Indeed, (3,3) maser emission has been detected 
previously in several regions of massive star formation mostly a ssociated with strong outflows 



e.g. PR 21(OH), NGC 63341, IRAS20126+4104. G5.89-0.39 Mangum and Woottenl Il994 



Kraemer and Jackson! Il995l ; IZhang et al.l Il999l ; iHunter et al.l 120081 ) . 



The stronger 25 GHz CH 3 OH maser spot (Fig. Eb) is coincident in positi on and velocity 



with one of the weaker 44 GHz Class I masers from iCyganowski et al.l (120091 ). and is within 
1.3" of the strongest 44 GHz maser. The weaker 25 GHz maser spot is located towards the 
SE edge of the CM1 UCHII region. Class I CH 3 OH masers are thought to be excited by 
collisions in shocks, with the shock liberating methanol from dust grains and thus providing 
the necessary high column density. ISobolev et al.l ( 120051 ) suggest that 25 GHz Class I masers 



require densities of ~ 10 cm and temperatures of 75- 100 K, hotter and den ser than 



that required for 44 GHz class I masers. As demonstrated by lVoronkov et al.l ( 120061 ) detailed 
modeling of the ratio of brightness temperatures among the 25 GHz masing transitions can 
help to pinpoint the physical conditions. We hope to carry out such analysis for the ensemble 
of 25 GHz masers observed across our survey in the future. 



5. Summary 

We present the first results from a 1.3 cm survey of MYSOs using the EVLA, focusing 
on the EGO source G35.03+0.35. The new EVLA correlator allows us to simultaneously 
observe the NH 3 (1,1) to (6,6) transitions, four CH 3 OH transions, two H recombination 
lines, and continuum. We find three major regions of dense NH3 gas, two (the NE and SW 
cores) being coincident with the EGO, and the third located in the adjacent IRDC. Using 
complementary SMA 1.3 mm and VLA 3.6 cm continuum data, along with temperatures 
derived from the NH 3 emission, we find warm gas temperatures in the range 20-220 K, 
core masses in the 20-130 M Q range, and average H 2 column densities of several xlO 22 
cm -2 . Together these data reveal a massive protocluster in an early stage of formation with 
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members representing different phases of MYSO evolution concurrently. This work highlights 
the potential diagnostic power that the full survey will provide toward understanding the 
formation of massive stars and protoclusters. 

Based on analysis using the CDMS, JPL, and Splatalogue spectroscopic databases, 
NASA's Astrophysics Data System Bibliographic Services, and the SIMBAD database. 

REFERENCES 

Argon, A. L., Reid, M. J., & Menten, K. M. 2000, ApJS, 129, 159 

Brogan, C. L., Chandler, C. J., Hunter, T. R., Shirley, Y. L., & Sarma, A. P. 2007, ApJ, 
660, L133 

Brogan, C. L., Hunter, T. R., Cyganowski, C. J., Indebetouw, R., Beuther, H., Menten, 
K. M., & Thorwirth, S. 2009, ApJ, 707, 1 

Cyganowski, C. J., Brogan, C. L., Hunter, T. R., and Churchwell, E.: 2011, ApJ 00 

Cyganowski, C. J., Brogan, C. L., Hunter, T. R., and Churchwell, E.: 2009, ApJ 702, 1615 

Cyganowski, C. J., et al. 2008, AJ 136, 2391 

Forster, J. R. and Caswell, J. L.: 1999, A&AS 137, 43 

Friesen, R. K., Di Francesco, J., Shirley, Y. L., & Myers, P. C. 2009, ApJ, 697, 1457 

Garay, C, Rodriguez, L. F., and van Gorkom, J. H.: 1986, ApJ 309, 553 

Hill, T., Burton, M. G., Minier, V., Thompson, M. A., Walsh, A. J., Hunt-Cunningham, M., 
& Garay, G. 2005, MNRAS, 363, 405 

Hunter, T. R., Brogan, C. L., Indebetouw, R., and Cyganowski, C. J.: 2008, ApJ 680, 1271 

Hunter, T. R., Brogan, C. L., Megeath, S. T., Menten, K. M., Beuther, H., & Thorwirth, S. 
2006, ApJ, 649, 888 

Kraemer, K. E. and Jackson, J. M.: 1995, ApJ 439, L9 

Kurtz, S., Churchwell, E., and Wood, D. O. S.: 1994, ApJS 91, 659 

Mangum, J. G. and Wootten, A.: 1994, ApJ 428, L33 



- 10 - 



Maret, S., Faure, A., Scifoni, E., & Wiesenfeld, L. 2009, MNRAS, 399, 425 

Perley, R.A., Chandler, C.J., Butler, B.J., Wrobel, J.M. 2011, ApJ, in press 

Rathborne, J. M., Jackson, J. M., and Simon, R.: 2006, ApJ 641, 389 

Rodon, J. A., Beuther, H., Megeath, S. T., & van der T ak, F. F. S. 2008, A&A, 490, 213 

Rosolowsky, E. W., Pineda, J. E., Foster, J. B., Borkin, M. A., Kauffmann, J., Caselli, P., 
Myers, P. C, & Goodman, A. A. 2008, ApJS, 175, 509 

Rosolowsky, E., et al. 2010, ApJS 188, 123 

Sobolev, A. M., Ostrovskii, A. B., Kirsanova, M. S., Shelemei, O. V., Voronkov, M. A., and 
Malyshev, A. V.: 2005, in R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley (ed.), 
Massive Star Birth: A Crossroads of Astrophysics, Vol. 227 of IAU Symposium, pp 
174-179 

Voronkov, M., Sobolev, A., Ellingsen, S., Ostrovskii, A., & Alakoz, A. 2005, Ap&SS, 295, 
217 

Voronkov, M. A., Brooks, K. J., Sobolev, A. M., Ellingsen, S. P., Ostrovskii, A. B., and 
Caswell, J. L.: 2006, MNRAS 373, 411 

Walmsley, C. M. and Ungerechts, H.: 1983, A&A 122, 164 

Zhang, Q., et al. 2007, ApJ, 658, 1152 

Zhang, Q., Hunter, T. R., Sridharan, T. K., & Cesaroni, R. 1999, ApJ, 527, L117 
Zinnecker, H. and Yorke, H. W.: 2007, ARA&A 45, 481 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 11 - 



Table 1. Observing parameters 



Parameter 



Value 



EVLA 1.3 cm (~ 24 GHz) AB1346 observations 



Observing date (duration) 
Configuration 
Primary beamsize 
Bandwidth 

Baseband subbands (MHz), a lines 



Baseband 1 subbands (MHz), a lines 



Velocity resolution 
Angular resolution 13 
Spectral line rms noise c 
Gain calibrator 

Bandpass and Flux Calibrator d 



07 Sep 2010 (3 hours) 

D 

2' 

16 x 8 MHz, single polarization 
23692.78, NH 3 (1,1) 
NH 3 (2,2) 



23724.78 
23828.78 
23868.78 
24084.78 
24140.78 
24508.78 
24532.78 
24927.88 
24959.88 
25023.88 
25055.88 
25295.88 
25327.88 
25687.88 



OH 2 n 9/2 (5- - 5+) 

NH 3 (3,3) 

S0 2 (8 2>6 - 9i, 9 ) 

NH 3 (4,4) 

H64a 

NH 3 (5,5) 

CH 3 OH (3 2) i - 3i, 2 ) 
CH 3 OH (5 2)3 - 5i, 4 ) 
NH 2 D (4 M - 4 0)4 ) 
NH 3 (6,6) 

CH 3 OH (8 2j6 — 8i j7 ) 
DC 3 N (3-2) 
H63a 

CH 3 OH (10 2)8 - 10i )9 ) 



25979.88 
0.4 km s -1 

3"7 x 3"0 (P. A. =-50°) 
3 mJy beam -1 channel -1 
J1815+0035 

J1924-2914 (17.1 Jy assumed) 
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Table 1 — Continued 



Parameter 


Value 




SMA 1.3 mm (~ 


225 GHz) observations 


Observing date (duration) 


24 Jun 2008 (11 hours) 


Configuration 


compact-north 


Primary beamsize 


52" 


Bandwidth 


2x2 GHz, single polarization 


Velocity resolution 


1.1 km s -1 


Angular resolution 6 


3"2 x 1"8 (PA.=+70°) 


Continuum rms noise 


3 mJy beam -1 


Gain Calibrators 


J1733-130 & J1751+096 


Bandpass Calibrator 


3C454.3 


Flux Calibrator 


Uranus 



a These are the rest frame subband center frequencies; be- 
cause the subbands were required to be on an 8 MHz grid, the 
lines are offset from the centers. Line rest frequencies were 
obtained from http:/ /splatalogue.net 

b Briggs weighting of 0.5. 

c Due to its inadvertent proximity to a subband filter edge, 
the NH 3 (2,2) subband is a factor of two noisier than the 
others. 

d 3C 286 observation failed, used bootstrapped flux density 
of J1924-2914 from previous track. 
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Positions for these fits are (J2000): CM1: 18 h 54 m 00.488 s , +02°01'18'.'21; CM2: 18 h 54 m 00.650 s , +02 o 01'19'.'05; NW CM2: 18 h 54 m 00.56F, +02°01'21('72; SW 
core: 18 h 54 m 00.164 s , +02°01'11'.'78; Dark Core: 18 h 53 m 59.396 s ,+02°01'01'.'47. 
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Fig. 1. — (a) Spitzer GLIMPSE three-color image of EGO G35.03 with RGB mapped to 
8.0, 4.5, and 3.6 fim overlaid with the EVLA primary beam (yellow dashed circle). Yel- 
low contours show the NH3 (1,1) integrated intensity (contours at 16, 28, 56, and 102 
mJy beam _1 *km s _1 ). The three brightest NH3 cores are labeled for reference. Magenta 
+ sy mbols and a single black o mark the locations of the 44 and 6.7 GHz CH 3 OH masers 
from ICyganowski et al.l (120091 ) . (b) and (c) show the NH 3 (1,1) velocity field (moment 1) 
and velocity dispersion (moment 2) with the lowest and highest (1,1) integrated intensity 
contour levels from (a) superposed in white. In all three panels, the EVLA synthesized beam 
is shown in the lower left. 
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Fig. 2. — Zoomed- in view of the NE and SW cores showing the NH 3 (1,1) integrated in- 
tensit y in colorscale wi t h VL A 3.6 cm continuum contours (0.1, 0.3, 2, 7 mJy beam -1 ) 
from ICyganowski et al.l ( 1201 ll ) superposed in blue; the three brightest centimeter contin- 
uum sources are labeled CM1, CM2, and CM3. (a) SMA 1.3 mm continuum is shown 
in white contours (levels: 7, 30, and 90 mJy beam -1 ) and the synthesized beams of the 
EVLA (yellow), VLA (blue), and SMA (white) are shown in the lower left corner, (b) 
White contours show the integrated intensity of NH 3 (2,2) (levels: 25, 50, 75, and 117 
mJy beam _1 *km s _1 ). (c) Contours show the integrated intensity of NH 3 (3,3) (white; 30, 
60, and 120 mJy beam _1 *km s _1 ) and (6,6) (cyan; 13, 17, and 22 mJy beam _1 *km s _1 ). In 
all panels, the locations of two new 25 GHz CH3OH masers are cyan triangles, and the 6.7 
and 44 GHz CH 3 OH masers from Fig [1] are also shown. 
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Fig. 3. — (a) Colorscale of the NH3 kinetic temperature derived from single Gaussian two- 
level fits; grey areas indicate where the assumptions for this technique are invalid. The blue, 
red, and white contours show the lowest 3.6 cm contour from Fig. [2j highest (3,3) contour 
from Fig. [2b, and highest (1,1) contour from Fig. QJi, respectively. The locations of the 
profiles shown in (b), for the non-thermal NH 3 emission (18 h 54 m 01.119 s , +02°01'14"08) and 
(c) the stronger 25 GHz CH 3 OH maser (18 h 54 m 01.076 s , +02°01 / 16"47) are indicated along 
with other positions fitted with method (2) and presented in Table El On (b) and (c), the 
vertical dotted line marks the LSR velocity of +53.1 km s _1 . (d)-(h) Observed NH3 spectra 
(black) overlaid with fitted profiles (red) using full LTE modeling in CASSIS (see Table [2]). 



